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Real-Space Renormalization Group for Spectral Properties of Hierarchical Networks 

Stefan Boettcher and Shanshan Li 
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We derive the determinant of the Laplacian for the Hanoi networks and use it to determine their 
number of spanning trees (or graph complexity) asymptotically. While spanning trees generally 
proliferate with increasing average degree, the results show that modifications within the basic 
patterns of design of these hierarchical networks can lead to significant variations in their complexity. 

To this end, we develop renormalization group methods to obtain recursion equations from which 
many spectral properties can be obtained. This provides the basis for future applications to explore 
the physics of several dynamic processes. 


I. INTRODUCTION 

Understanding the dynamics of hierarchical systems is 
rapidly becoming a subject of wide interest within the 
study of complex networks This may seem surpris¬ 

ing, as hierarchies are not a new subject in complexity [^, 
and considering that the ultimate hierarchical system - 
a tree - has long been studied in many variations and 
for many types of dynamics because exact results can 
be obtained. Known as Bethe approximation in statis¬ 
tical physics BIS], trees locally resemble sparse mean- 
field models in physics, due to their infinite dimension¬ 
ality. Beyond the mean-field statistical models, however, 
there exists a much richer dynamics arises from peculiar 
structural features of hierarchical networks such as high 
degree of clustering and modularity. In many realistic 
situations, such hierarchies are embedded in some finite¬ 
dimensional space mio], which can lead to an entirely 
new set of phenomena. Examples of such embedded hier¬ 
archical systems contain transport and control systems, 
they have been observed in the brain mm. or they 
have been studied for their novel synthetic critical behav¬ 
ior [IMT], such as explosive transitions that are purely 
induced by the geometry in percolation m or synchro¬ 
nization m- Recently, hierarchical networks based on 
Dyson’s model have shown to allow for the existence of a 
number of metastable states, which can be used to study 
the modular architecture and parallel processing in neu¬ 
ron networks as well as the ergodicity breakdown for the 
stochastic process 

Key to understanding the mechanisms at the core 
of those novel phenomena lies within the geometry of 
these hierarchies. These geometric properties are in¬ 
evitably tied to spectral properties of their network 
Laplacians m- Here, we investigate a recently proposed 
class of hierarchical networks [1] within the simplest of 
spatial embeddings - a simple line - for which we can 
obtain many spectral properties exactly using the renor¬ 
malization group (RG) [7] .For these Hanoi networks, we 
study some basic properties of their Laplacians. In par¬ 
ticular, we derive sets of recursion relations that allow to 
study their secular equations (also called characteristic 
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polynomials), whose zeros provide all eigenvalues, to ar¬ 
bitrary accuracy. We determine the asymptotic scaling of 
their determinants with system size. These determinants 
can be used as a generator for many other asymptotic 
properties. 

The Laplacian matrix is given by 

= diSij — Aij, ( 1 ) 

where di specifies the degree of the /-th site and Aij is 
the adjacency matrix of the network. Since we assume 
that the links in the networks are undirected, A, and 
hence L, are symmetric. We further assume that there 
are no external links, which implies the vanishing of all 

row or column sums in L, The 

fundamental property characterizing the Laplacian ma¬ 
trix is its spectrum of eigenvalues, the solutions X; of the 
secular equation 

det[L-Al] = 0. (2) 

This spectrum is highly nontrivial for the Hanoi net¬ 
works, and we will not be able to describe it in much 
detail here. But we can provide an RG approach that 
reduces the effort exponentially from solving 2^ x 2^ de¬ 
terminants to k iterations of a closed set of RG recursion 
equations for any desired quantity, where N = 2^ refers 
to the number of sites in the network. Numerous aspects 
can be extracted in closed form asymptotically. 

The spectrum of network Laplacians features in many 
practical applications. The scaling of the ratio between 
largest and smallest eigenvalue indicates the synchroniz- 
ability of a network m. which also can be approximated 
by a sum over all eigenvalues (inverted) when related to 
a random deposition process [2^. Permeability and 
well-connectedness of networks can be defined in terms of 
their smallest eigenvalue m- The spectrum further fea¬ 
tures prominently in quantum transport m ,the behavior 
of continuous-time quantum search algorithms [291 ld(J] , 
graph partitioning mm and image processing [55], iust 
to name a few examples. Of course, Laplacian spectra 
also determine the characteristic frequencies of mechan¬ 
ical vibrations, from which connectivity between inter¬ 
acting units within the membrane can be identified |34j . 
Thus, there has long been strong motivation to study 
such spectra particularly on fractal networks, where its 
properties can be explored in great detail [3SJ3D]. 
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Figure 1. Depiction of the Hanoi networks HNS and HNS. 
The 3-regular network HNS corresponds to the solid (black) 
lines alone, while HNS in addition also consists of the (green-) 
shaded lines. For HNS, sites on the lowest level of the hierar¬ 
chy have degree 3, then degree S, 7, etc, concerning a fraction 
of 1/2, 1/4, 1/8, etc., of all sites, which makes for an average 
degree S in this network. Note that both networks are planar. 



Figure 2. Depiction of the nonplanar Hanoi networks HNNP 
and HN6. Again, starting from a Id-backbone (black lines), in 
HNNP a set of long-range links (blues-shades lines) is added 
that break planarity but maintain the hierarchical pattern set 
out in Eq. ([^ that leads to a network of average degree 4. If 
we add again the same links that distinguished HNS and HNS 
in Fig. we obtain HN6 with average degree 6. In all cases, 
the RG on these networks remains exact. 


Here, we will only focus on the simplest case of the 
scaling of the determinant itself, obtained by taking A —>■ 
0 in Eq. ([^ , which provides the number of spanning trees, 


#ST = - lim 
A->.0 


det [L — Al] 
XN ^ 


( 3 ) 


sometimes also referred to as “graph complexity ” mi- 
Eq. ^ is one of the oldest results in algebraic graph the¬ 
ory, due to Kirchhoff (1847) [ 13 ] • Spanning trees describe 
the size of the attractor state in the self-organized critical 
sandpile model jUj, they characterize the optimal paths 
between any two sites in a network |43j . and are also re¬ 
lated to optimal synchronizability of a network |44] . The 
number of spanning trees is of fundamental interest in 


mathematics and physics. For example, it is related to 
the partition function of ( 7 -state Potts model in the limit 
q^O gUSS] Thus, studies on the asymptotic growth 
of spanning trees are well motivated not only for regular 
lattices m, but also on self-similar structures gHHSOj. 
On these networks, the number of spanning trees =ffsT 
exponentially increases with N, which can be character¬ 
ized by the tree entropy, which is the entropy-density of 
spanning trees miiaisi], 


s = 


lim 

N^OO 


ln(#ST) 

N 


( 4 ) 


The number of spanning trees on Hanoi networks can 
be easily derived from our RG procedure, allowing us to 
explore the role of geometric structures on its asymptotic 
growth. 


II. STRUCTURE OF HIERARCHICAL 
NETWORKS 

The Hanoi networks we discuss in this paper were in¬ 
troduced first in Ref. gj. Each of the networks consid¬ 
ered in this paper possesses a simple geometric backbone, 
a one-dimensional line of sites N. in which each site is 
at least connected to its nearest neighbor left and right 
on the backbone. To generate the small-world hierarchy 
in these networks, consider parameterizing any number 
n> 0 uniquely in terms of two other integers (*,/), i ^ 1 
and 1 < j < 2 ^“*, via 

n = 2*-i(2j-l), (5) 

which motivated the name of the networks [52|. 

The networks can be considered either closed into a 
loop of = 2 ^ sites or in a linear arrangement with 
N = 2^ + 1 sites, a difference that typically does not effect 
the RG and any leading asymptotic property. Their re¬ 
cursive pattern is far more easily to illustrate when drawn 
linearly, as for HNS and HN5 in Fig. [^and for HNNP 
and HN 6 in Fig. although in our calculations we avoid 
spurious edge-effects by considering periodic loops, where 
sites n = 0 and n = N = 2^ are identical to each other. 
For convenient comparisons, we call the ordinary one¬ 
dimensional loop HN2 (for Hanoi Network of degree 2). 
The details of the design and most of their geometric 
properties have been discussed at length elsewhere m- 


HI. OBTAINING THE RG-RECURSIONS 

In this Section, we show how to obtain the recursion 
equations for the RG-flow of our networks that describes 
the asymptotic properties of the lattice Laplacian. To 
that end, we employ the well-known formal identity |5S] , 
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1 

\/det4“^ 


CXD 



— CxD 



( 6 ) 


for the evaluation of the determinant of for each net¬ 
work HNa with a G {2,3,5, 6} of size N = 2^. (Since 

Laplacians are singular matrices, formally, some integrals 
will not converge and need regulation, however, these sin¬ 
gularities will be dealt with explicitly below in the final 
step of the RG.) We can group our networks into two 
classes of equivalent topologies, as respectively Figs. 
and [^suggest: HN3 is contained within HNS and HNNP 
is contained within HN6, while the simple loop (HN2) is 
contained within both, of course. In addition, removing 
the backbone of HNNP decomposes it into three sepa¬ 
rate, loop-less tree structures that have been studied in 

Ref. [51] . 

The following calculations are purely formal, as some 
of the integrals may not converge in Eq. ([^, in which 
case we would have to insert some form of regulation. 
Yet, the integrals merely serve as generators for the de¬ 
sired RG-recursions, essentially. Alternative algebraic 
means to evaluate the determinants that do not involve 
integrals have been used for a specific application pre¬ 
viously |1S1 [35]. However, the approach taken here 
provides entirely equivalent results with a clearer per¬ 
spective on the topological transformations that are in¬ 
volved. Furthermore, the current approach is closely 
linked to the description of a path-integral for a free 
field-theory ESI on a given network, which affords cer¬ 
tain extensions of the methods, for example, by inserting 


source-terms into the exponential in Eq. ^ 

as generators for more complicated observables. When 
these sources are either uniform {Jn = J f.a. n), local¬ 
ized {Jn = JSn,x), or hierarchically staggered [Jn = Ji(n) 
according to Eq. §], exact RG-recursions may still be 
obtained. 


A. Renormalization Group calculation for the 
determinant of HNS and HNS 

Here, we evaluate the most general 2^ x 2^- 
determinant for the matrix that is used in the anal¬ 
ysis of HN3 and HN5 (and HN2) below. The prop¬ 
erties of each determinant emerges via renormalization 
from for a different set of “bare” parameters, while 
the RG-recursions themselves remain the same for each 
a G {2,3,5}. These parameters describing the renormal¬ 
ized weights of effective links between sites, while their 
bare values serve as initial conditions on these RG recur- 
sion equations. We define L} ' as the matrix of a one¬ 
dimensional loop of 2^ sites (which we may call HN2). 

and are respectively the matrices for the HN3 
and HN5 networks. Instead of providing a formal de¬ 
scription of for general size 2^, we simply illustrate 
its generic recursive pattern for the case fc = 4: 
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The bare parameters qt on the diagonal refer to on-site 
properties of each site n that belongs to the Fth hierar¬ 
chy as determined by Eq. ([^ . For the case of the lattice 


Laplacian here, for example, qi is simply the degree of 
that site. The off-diagonal parameters label extant or 
potentially emerging links between sites that are undi- 
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rected, so the matrix is symmetric. To keep parameters 
fundamentally non-negative, we insert negative signs ex¬ 
plicitly, so that the bare values for all pi is uniformly 
unity, however, they renormalize differently for each i. 
(Here, level * = 0 refers to nearest-neighbor links in the 
backbone.) While the pi correspond to the solid black 
lines in Fig.[2 the parameters li refer to those long-range 
links shaded in green in Fig. whose addition make up 
HNS. (Entries like —21^-1 near the highest level of the 
hierarchy correspond to a convenient choice in imposing 


periodic boundary conditions on the network.) Although 
all li are originally zero for HNS, we still need to consider 
them, as they emerge as a relevant parameter [ 7 ] during 
the RG, even for HNS. If p* = 0 for all i > 0, Eq. Q 
reduces to the tridiagonal matrix for a simple loop we 
call HN2. Note that we have imposed periodic bound¬ 
ary conditions by identifying site N = 2^ with site n = 0, 
where the site indices run from n = 0ton = 2^ — 1. 

Employing the binary decomposition of the integer la¬ 
bels implied by Eq. ([^ for the sites on the network back¬ 
bone, we can write Eq. as 


detL 


(3,5) 


■ 2 _ ^ JJ dxodx2k- 

— OO 

CxD ^ fc— 1 

" nn 

i=i j=i 


-exp{-qfe (xl + xlk-i) + 4:lk-iXoX2k-i} 


( 8 ) 


dx2i-l(^2j-l) 


fe -1 2 '=“* fe -2 2 '““‘ 

~ ^2»-l(2j-l) ^2»-l(2j-2)2^2»-l(2j) 

i=l j = l i=l J = 1 

k-l 2'“-! 

+2 E Pi ^ —3)^2^“l(4j —1) ^ ^ ^2jr —1 (^2^ —2 ^2j) 

i=l j = l j = l 


The factor 1^, initially unity, captures the contribution 
of integrals from any prior RG-step. 

(3 5) 

To solve deti), ’ ' recursively, we integrate only 
over all variables x of odd index [those with i = 
1 in Eq. ([^]. To that end, we focus on the 


case i = 1 in the product of integrals in Eq. 

b 2 i-i) 


_ 1 - 2 ^ _~ 

and re-write 117=1 '^^2»-i(27-i) , == 11^=1 dx2j-i = 


1=1 


11^=1 dx 4 j- 3 dx 4 j-i to get: 


) K — ^ LX^ 

JJ" // ——MlLLexp{-gi (a;4j_3-|-a;4^_i)-l-2pia;4j-3a;4j-i -l-2po [a;4i-3 (a;4j-4+ a:4j-2)+ a;4j-i (a:4j-2+ 2 : 4 ^)]}, 


j=i- 


= n (9?-P?)“"exp 

j=i 


giPo 

qf-pl L 


= {ql-pl)-^ exp 


{X4j-4: +X4j-2) + {x4j-2 + X4j) +^2-^ (2^4j-4 + 2:44-2) (2:44-2 + a:4j) 

qi~Pi 


(9) 


2 giPo I 2 1 I ^qiPo 

2 _ 2 V-^ 0 -l“ 2 , 2 fc_ij-I- 2 _ 2 ^ ^ 

7=2 7 = 1 


E E ^ 2 »- 1 ( 24 - 1)2 


i-)k — 2 

2 Po V -.„2 


gi-pi ^ qi-pi . 

1=1 j=i 


I ‘^Po I 2piPg i 

+ „ 2:24-22:24 + -2 y2 }_^ X(^2j-2)2X(2j)2 ^ • 

gl Pi 4 = 1 


With that result, the remaining integral over the even- 
indexed variables can be written as 
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detL 


(3.5) 


II ?'E2!^exp 


.2 1 


2 giPo 

yfc 2 2 

9!-PI 


(a;o + a::2''-i) +4^fe-i^oa;2fc-i 


/A; — 2 2 ^ 

n n 

-oo \*=1 2 = 1 


da; 


2 d 2 j-l) 

2 pg 


fe -2 


eiiPS -E 


i=2 


Qi-\-l ■ 


2 giPo 
^2 ^2 


E 


'' 2 »( 2 jr-l) 


92- 

91-P1 
fc-3 2 


2““^ k—2 2'““*“^ 

E ^2(2j-1)+2Ep*+ 1 E ^2i(4i-3)2^2*(4j-l) 


( 10 ) 


4 = 1 
fc —i — 1 


j=i j=i 


2 \ 2 

+2 E h+l E ^2»(2i-2)^2d2i~)+2 I ^2+ 2 ^2 ) E ^2(2j-2)^2(27) 


i=2 j = l 




+21 + ■ 


2 \ 2 
Po 


91-Pi 


^ E ^2(2j-l) (2;2(2i-2) +a;2(2j)) 


Substituting x'^ = a; 2 i, this expression is identical in form 
with Eq. §, and we can identify 


9i =92-2- 


Po 


91 -Pi 


/ _ . o ' 2 'iPo 
Qi Qi-\-l ^2 2 ’ 

Ql-Pi 


p'o — h + ■ 


Po 


91 -Pi 




Pi =Pi-hi> 

7 / _, , PiPo 

11 ^2 “T 2 2 ’ 

91-PI 

^i = ^j-l-i, (*> 2 ), 


(*> 2 ), 


( 11 ) 


The difference between the primed and unprimed quan¬ 
tities represents the step from the /i-th to the /a -I- 1-st 


level in the RG recursion, with 0 < ^ < k. The recursion 
for the overall scale-factor requires more care, as it 
depends explicitly on k and we have to take into account 
the level /a at which the factor in front of the integral in 
Eq. (10) arises, thereby shifting k^ k — fi. Thus, 


I'k = lk{ql-Pl) 


_ nk — fi—3 


( 12 ) 


B. Renormalization Group calculation for the 
determinant of HNNP and HN6 


Now, we evaluate the most general 2^ x 2^-determinant 
for the matrix representing either the HNNP or 

the HN6 network. Again, we simply provide a description 
of for the case fc = 4: 
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(13) 


The meaning of the bare form of the renormalizing pa¬ 
rameters (qi^pi,li) is the same as in Sec. Ill A In par¬ 


ticular, the dark-shaded links in Fig. [^correspond to pi, 




























6 


while the green-shaded ones again refer to h, which may in Eq. (|^ to write 
be originally absent but emergent in HNNP. The deter¬ 
minant of can be evaluated by using the identity 
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Substituting back into Eq. (141 obtains: 
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Substituting x'^ = X 2 i, this expression is identical in form 
with Eq. (141, and we can identify 


qi 
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I' = I 


The difference between the primed and unprimed quanti¬ 
ties represents the step from the ^-th to the p, -I- 1-st level 
in the RG recursion, with 0 < ^ < k. As for HNS above, 
the recursion for the overall scale-factor requires a shift 
k ^ k — 


Ik — Ik 9i 


(18) 


where we used that and = 0 for all 


i > 1 while Pfi = 


The recursion for in either of 


Eqs. (12) or (18 1 has the formal solution 


An) 


n 

i=0 


ry n, 

9r 


( 21 ) 


After the k — 1-fold application of the RG-recursions in 
Eq. ( |20| reduces the original 2^ x 2^ matrix - such as in 
Eq. - down to a 2 x 2 matrix for a “loop” with only 
two sites that are doubly linked, and we have 


det 


Lf-Al 


,(fe 1) 


-2 


det 


9fc-i 

-‘IPk-i 


‘Ipk—l^ 

9fe-i 


2 ) 


Exact Solution for HN2: The nonlinear recursions in 
Eq. ( [^ are easily solved in closed form by defining = 
q 11 /Pill ^°r which — 2 , obtained by dividing the 

2nd by the 3rd line. The solution is 


s^ = 2 cos 


2 '^arccos 

2 poy 


= 2 T..(|^) (23) 


IV. 


RG-EVALUATION OF NETWORK 
LAPLACIANS 


Here, we use the RG-recursions in Eqs. (11 1 and rtl7| to 


determine the asymptotic scaling behavior of the deter¬ 
minant of the network Laplacians for large system sizes. 
We begin with the example of a simple line, HN2, which 
is contained in either equation, to re-derive some familiar 
results to demonstrate the procedure. 


A. Simple Example: One-dimensional Lattices 

The RG allows us to find the secular equation for HN2, 
a one-dimensional loop oi N = 2^ sites, as a reference. In 
that case, all sites have constant degree di = 2. With 
the RG approach from Sec. M we have to solve the 
recursions either in Eqs. ( TT|l 2 ) or in Eqs. ( TTflS l for 
the initial conditions on the bare parameters. 


4°^ = 1 

gf = 2-A 
=1 
P?^ = 0 


(i^i) 


(19) 


= 0 


(i^l) 


Note that we allowed here for a prospective eigenvalue A 
subtracted from each diagonal element, as indicated by 
the eigenvalue Eq. (i). In both sets of RG-recursions, the 
equations simplify to 


9m+i = 9m 




9m 


Pm +1 = 


9m 


( 20 ) 


where Tn{x) refers to the n-th Ghebyshev polynomial of 
the first kind [SS]. Inserting into Eqs. (20 1 and applying 


the initial conditions in Eqs. (191 generates the results 

Pm = n 77’ 9 m = s 

j=0 


^M 


M-l 

n 


1 


M-l 


1 


i=0 


= n 7 

• n 
2=0 


where the last equality emerges from Eq. (21 1 under 


reordering factors in the products. Alternatively, one 
could realize that the Ansatz 


Pm = 


sC 


2 M 


1 -C 


2M + 1 ' 


9m = 




i-C 


2M+1 


(25) 


also provides an exact solution of Eq. (20 1 , which match 


the initial conditions at /x = 0 for s = Ei ^(4 — A) A and 
C = 1 — ^ =F 2 \/( 4 — A) A. In either case, note that for 
A —>■ 0 the solution reduces to Ifi=Pfi = 9m/2 = 2~^. 

We note that these formal solutions, albeit closed-form, 
are rather complicated and even numerically very difficult 
to us e fo r arbitrary A. However, inserting Eq. (241 into 
Eq. (22 1 with N = 2^ provides 


det 


l(^)-ai 


= 2T^(1-- )-2, 


(26) 


a well-known exact result. Expanding Eq. (26 1 to first 


order in A provides for the number of spanning “trees” in 
Eq. (§: 


A2)_ r^(l) 

#ST- ^ 


= iV, 


(27) 


which is exactly the number of open strings covering all 
sites that one can embed on a closed loop of N sites. 

































Asymptotic Solution for HN2: In general, we will not 
be able to solve the nonlinear recursion equations of the 
RG-flow in closed form, of course. It is therefore instruc¬ 
tive to explore asymptotic methods for such cases by way 
of this exactly solvable instance. We note two proper¬ 
ties of the RG-flow that turn out to be quite general: 
(1) Evolving numerically from the initial, bare parameter 
values in Eq. (19) for A = 0, these parameters approach 
their asymptotic value at ^ > oo with a correction that 
decays as for some a > 1 ; and ( 2 ), to remove the triv¬ 
ial A = 0 eigenvalue that each of the network Laplacians 
Lk possesses, we further need to expand the RG-flow re¬ 
cursions to first order in A for A —0 while /i —>■ oo . As the 
stable fixed point for /i —>■ oo at A = 0 becomes unstable 
for A > 0, that correction diverges with some power /3^ 
for /3 > 1. Therefore, we need to make an Ansatz typical 
to analyze unstable critical points in RG mm- 

Pii"^Poo + oi ^Po +XI3^Pi +..., 

Q^i'^Qoo + ct ^Qo +XP^Qi +..., (28) 

for /i —> oo and A —> 0. For HN2, poo = Qoo = 0. Then, we 
obtain at A = 0 and to leading order in a~^: 


a Qo 

which have the solution 


Po 


p2 

fJL 

Qo’ 


a — 2, Qq — 2Pq. 

Extending to the order-A correction yields 


PQi = -Qi-2Pi, 


with solution 


^ = 2 , 


PPi=Pi-\Qi, 


Qi = -4Pi. 


(29) 


(30) 


(31) 


(32) 


The first line in Eq. (20 1 in that form is difficult 
to evaluate. An asymptotic evaluation doesn’t seem 
possible, generally, since most contributions arise from 
the terms with smallest /r where the power 2 ^“^ is 
largest but in Eq. (28 1 has not achieved its asymp¬ 
totic form yet. It is easy to evaluate numerically to 
any accuracy, though, when rewritten as a recursion in 

k, • Here, one easily finds that 

= 2 “* = 1 /A^, since = Pk-i for all k, of 

course. Here, no 0(A)-correction was needed, since the 
denominator of Eq. (22 1 is already finite at A = 0. In con¬ 


trast, the numerator cancels at A = 0 and yields instead 
q‘l_i —^Pk-i ~ ^ (q) to next order. We Anally get 


for Eq. (22 I: 


det 


Lf-Al 




(33) 


where we have ignored overall factors. In this simple 
case, it happens to be a = (3. Taking these facts into ac¬ 
count, Eq. (33) reproduces the exact result in Eq. (27). 


Note that to obtain the dominant contribution (if it ex¬ 
ists) that varies exponentially in system size N = 2^ 


we 

merely need for A = 0; the remaining integral con¬ 
tributes at most a power-law pre-factor, aside from the 
A = 0 trivial eigenvalue itself. 


B. Case HNS 

In this case, we have to interpret the results in 
Eqs. © for the bare parameters 


j(o) _ ^ 

— -Li 


9* 


( 0 ) _ 


= 3-A 


( 0 ) 1 
Pi’ = l 

l?^ = 0 


(*^ 1 ), 

(*^ 0 ), 

(*^ 1 ), 


(34) 


reflecting the fact that all sites in HN3 have a constant 
degree di = 3. Since all diagonal entries are identical, the 
hierarchy for the qi collapses and we retain only two non¬ 
trivial relations, one for qi and one for all other qi = q 2 
for all j > 2. Here, all pi are non-zero, encompassing the 
backbone links {i = 0 ) and all levels of long-range links 
(j > 1). But it remains pi = 1 for i > 1 at any step p, 
of the RG, in particular, = 1 throughout; only the 
backbone po renormalizes nontrivially. Although all bare 
links of type li are absent in this network, the details of 
the calculation in Sec. IHI Al show that under renormal¬ 
ization terms of type li emerge while those for li for i>2 
remain zero at any step. Thus, Eqs. ( II|l 2 ) reduce to 
RG recursion equations that are far more elaborate than 
for HN2 in Eq. (20) above: 


7 -( m + i ) _ An) 
— -'fc 


in) 

q\ 


di 


(m+ 1 ) _ An) 9 L 
1 — ^2 ^ 


in) 

Po 


-1 

2 






q(/+P = qM-2- 


An) 

r (/a )1 

2 

9i 

[Po 



in) 

q\ 


2 ’ 
-1 


(35) 


Po 


(m+ 1 ) _ lin) I 
0 — A + , 


An) 

Po 




‘1 


in) 

Po 




We can further eliminate q 2 from Eqs. (35) by noting 
that these equations possess an invariant: 




-21); 


in) 


{0 A p <k). 


(36) 
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Then, abbreviating = q[^\ = p\^\ and 

Eqs. (351 reduce to 


Q/j+i — Q'/j “t" 2/^ 2 


pI 


(go = 3-A), 


P^J,+l — + 


9m 


-1 


{Po = 1), 


(37) 


Pu 


7 _ ^ M 

9^-1 


(^0 = 0 ). 


We can again solve for in Eq. (351 independently, 

w-i] 


1 —2 


(38) 


j=0 


Evolving the RG for k — 2 steps results in a reduced 
network that consists of a loop of 4 sites, formerly at 0, 
2^-^, 2*“^, and 32^“^. The sites at 0 and at are 
now connected doubly by links lk -2 whereas the other 
two are connected by a previously unrenormalized link of 
bare unit weight. Each site is of course connected to its 
nearest neighbor in the backbone loop with renormalized 
links Pk- 2 - Note, that the sites at 0 and at 2^~^ have a 

on-site factor of q^^ = q[^ - 1 - 24 ^ = qk -2 + ‘^h -2 

whereas the other two sites have a factor of q[^ = qk- 2 - 

Hence, the secular determinant for HN3 reads: 


det 




j(k 2 ) 


r(fe- 2 ) 


^ 2 


det 


9fe-2 + 24-2 —Pk-2 —2Za;-2 —1 —Pk-2 

—Pk-2 Qk-2 —Pk-2 —1 

— 24-2 — 1 —Pk-2 9fe-2 + 24-2 —Pk-2 

—Pk-2 —1 —Pk-2 qk-2 


-2 


{qk -2 + 1 ) {qk -2 + 44-2 + 1 ) \{qk -2 — 1 )^ — 4 Pfe_ 2 ] 


(39) 


The evaluation of Eq. (391 follows closely the analysis 
in Sec. |IV A[ Our numerical investigations indicate that 
the recursion equations in (37), starting from the initial 
conditions with A = 0, evolve to a fixed point at which 
goo = 1 while Pfj, Ifi ^ Poo = loo = 0. Thus, an 


Ansatz for fixed points |7] similar to Eq. (28), 


g^"^l + a ^Qo + A/3^(5i, 
Pfi ~ a“^Po +A/3^Pi, 
4^0! ^Lo + A/3^Ti, 


(40) 


for /i —>■ oo and A —>■ 0 and requiring a, (3 > 1, when in¬ 


serted in Eq. (37), yields the unique solutions 


_ 2 Qo 

“ (i’ 2 ’ 


/3 = 2 , Pi = - 


5Qi 

12 


j Qo 

, Li = -^. 

6 


(41) 


Here, (j)= (■\/5 + l) /2 = 1.618... is the “golden ratio” [57] . 
and Qo,i remain as arbitrary overall constants, whose 
knowledge would require a global solution of Eqs. (37). 
When applied to Eq. 

becomes ^ ^ XN. 


nant 


(39), 


becomes 


Ak-2) 


-2 


the determi- 
The pre-factor 


-\k — 2 — II 


= n:i^oK-i] 

-f+2(fe-l) 


yk — 2 — IJ, 


n*=o ~ + , where x 

is a constant that can be determined to any accuracy 
from Eq. (|38|). For instance, it provides the recursion 

XT 2 


j{k 1 ) 

^k + 1 


,-l 


r(fe-2) 


extract x 


r{k—2) 


n- 2 'l 1 /^ 


Jk-2) 


-2^ 


for 


large k, e.g., x « 2.0189990298 at fc = 40. Inserted into 
Eq. (39), we obtain 


det 


4 ")-ai 


' AAf2-log2</,^Ar 


(42) 


for A —>■ 0, where we have ignored any pre-factors again. 
Note the similarity to the calculation leading to Eq. (26). 
ByEq.(§, we then find for the number of spanning trees 
on HN3: 


n(3) 

TTST 


^ 1-1082 02 , 0189990298 ^. 


(43) 


HNS without backbone: An extreme check for the con¬ 
sistency of the RG recursions in Eqs. ( II|l2 ) is provided 
by the degenerate case of Hanoi networks without back¬ 
bone. For HN3, t his w ould imply that the bare param¬ 


eter equations in ( |34[ ) are modified to g-°^ = 1 — A and 
Po^^ — case, all recursions in Eqs. ( ll][l 8 ) be¬ 

come trivial, with pjf'^ = l^^^ = 0 and q^^^ = 1 — A for all 
z > 0 and 0 < p < k while p[^'^ = 1 only for z > 1. Now, 


since qf—pf^ 2X for all p, we find If .^2 (—2A) ^ ^ and 
that the determinant of the last four sites merely has sin- 

N 

gle line with ~ —2A, such that det[Lk —Al] ~ (—2A)^. 
This correctly reflects the fact that HN3 without back¬ 
bone decomposes into ^ disconnected individual links, 
see Fig. 12 each link by itself has a trivial determinant 

^ ^ —2A. Glearly, each such determinant 


for k > 2 from which we = 


-1 1-A 
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divide by the number of its sites (= 2) and A, according 
to Eq. (|^, merely says that single links have a unique 
spanning tree. 


C. Case HNS 

HN5 contains HNS but possesses many additional links 
between sites, see Fig. [T] and therefore, we expect that 
the number of possible spanning trees proliferates faster 
in HN5. In this case, we have to interpret the results in 
Eq. (11) for the “bare” parameters as 


iW _ 1 

q^°'> = 2i+l-X (i^l), 

pP = 1 (i^O), 

Zf) = l (z^l), 


(44) 


reflecting the fact that all sites in HN5 have a hierarchy- 
dependent, increasing degree of d* = 2 j -|- 1 with average 
5. Now, diagonal entries are no longer identical and we 
have to modify the equations for the qi when compared 
to HNS. Yet, the difference between qi and is con¬ 
stant throughout, — qi = 2, and taking that modifi¬ 
cation into account, only the renormalization of qi and 
q 2 evolves nontrivially, as before for HNS. Again, all pi 
are non-zero, encompassing the backbone links (z = 0 ) 
and all levels of long-range links (i > 1). But it remains 
Pi = 1 for i > 1 at any step p of the RG-flow, in particular, 
= 1 throughout; only the backbone po renormalizes 
nontrivially. Special to HN5, all links of type li are al¬ 
ready present initially in this network. Though, only li 
renormalizes, as in HNS, whereas it is L = 1 for all i> 2. 


Thus, we obtain more elaborate RG recursion equations 


which merely differ in the last relation from Eq. (S5): 




= l-k- 


ik) 

Po 


1 2 


(m) 


(45) 


-1 


The formal solution for is unchanged from HNS, 
given in Eq. (S 8 ). Furthermore, we note that, despite 


the changes to the recursions for q 2 and /i, the invari¬ 
ant in Eq. (S 6 ) remains valid, allowing the elimination 

for p < k — 2. Then, abbreviat- 
reduces the 


An) 


of the recursion for qif 

in) in) 

mg again q^ = qf', p^ = p^f', 

RG-recursions to 


and = 4 ^^ 


Qfi+i — Q^i T 2/^ 2 


Pr 


<ln 


-1 


(go = 3-A), 


Pn+i 


^+1 


= 1 ,. 


= l-k 


<ln 


-1 


(po = 1 ), 


(46) 


pI 

9^-1 


(?o = l)- 


As in HNS, evolving the RG for k — 2 steps also results 
in a reduced network that consists of a loop of 4 sites, 
formerly at 0, 2*“^, 2^“^, and S2^“^. The sites at 0 
and at 2^~^ are now connected doubly by links lk- 2 - 
Each site is of course connected to its nearest neighbor in 
the backbone loop with renormalized links Pk- 2 - Note, 
since the invariant in Eq. (S 6 ) is invalid for p = k — 2 
due to degree dk^2.k + l for sites 0 and 2 *“^, special 
consideration is required for the on-site factor of q^ 
abbreviated as rfe_ 2 , whereas the other two sites have a 

factor of q^^ = qk- 2 - Hence, the secular determinant 

for HN5 reads: 


det 


l(")-ai 


jik 2 ) 


jik 2) 


det 


1 -2 


—Pk-2 —‘^h-2 —Pk-2 
—Pk-2 qk-2 —Pk-2 —1 

— 24-2 —Pk-2 rk-2 —Pk-2 
—Pk-2 —1 —Pk-2 qk-2 J 


(47) 


{qk -2 + 1 ) ( 24-2 + 'rk- 2 ) [(1 — qk- 2 ) ( 24-2 — 'rk- 2 ) — 4 pfe_ 2 ] 


The seemingly innocuous difference in the last relation 
of Eqs. (46) compared to Eq. ( |^ for HNS has dramatic 
consequences. Instead of the singular scaling Ansatz in 
Eq. (40) we used for HNS (or HN2), Eq. ( [4^ has an 
ordinary fixed point at /x —> 00 with goo = (5 + -\/4T)/2, 
Poo = 2/00 = (3-I- 4 ^ 1)74 as algebraic solutions at the 
fixed point. A simple perturbation for small A on Eq. 


I 

then yields [7] 
9 m" 
Pm" 
4- 


S-kVil 

2 

3-kV4T 

4 

3-kV4r 


-l- A2^Qi, 

-A2^Qi 

-A2^Qi 


57 - 7^41 

40 ’ 

47-7^44 
40 ■ 


(48) 


By the same method as described in Sec. IV B we ob- 

_ '2 

tain a; « 2.7548806715 at fc = 40 for 
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Inserted into Eq. (471, we obtain (up to a factor): 

AiV2.7548806715^, 


det 




or, for the number of spanning trees, 


rrsT 


2.7548806715 


N 


(49) 


(50) 


HNS without backbone: As Fig.j^suggests, HN5 with¬ 
out its backbone implies = 0 but Z= 1 and 

= 2i — 1 — A for all i > 1. Then, the first itera¬ 
tion of the RG recursions in Eqs. ( Tl|T 2 ) evolves triv¬ 
ially, producing p^'^ = 1 and = 2z -|- 1 — A 

as well as = 1 for all i > 1. If we also en¬ 

force = 1 , these are just the initial conditions of 


the RG-recursions in Eq. (441 again, except starting at 


p= I, and we reproduce the same result as in Eqs. (49 


501 but for a network of size N/2. One complication 
arises from the N/A disconnected pi-links. These are 


accounted for via the first iteration of in Eq. (|12|: 


with = 1 — A and p^^' = 1, it produces a factor of 


JO) 


= [qi—Pi] ~ (— 2 A)^, i.e., exactly one fac¬ 
tor of —2A for each disconnected line. In the calculation 
of the spanning trees of the remaining network, these 
need to be ignored, hence, requiring = 1 . 


D. Case HNNP 


2 < k, as for HN5 above, but instead with average 

degree 4. Here, too, the difference between qi and qi+i 
for * > 2 is constant throughout, qi^i — qi = 2, such that 
only the renormalization of qi and q 2 evolve nontrivially. 
Again, it remains pi = l for i ^ 1 at any step p of the RG, 
in particular, = 1 throughout and only the backbone 
po renormalizes nontrivially. Although all links of type 
li are initially absent in this network, under renormaliza¬ 
tion terms of type li emerge while those for li for * ^ 2 
remain zero at any step. These considerations reduce 
Eqs. ( [iTp^ to 


91 -92 ^ (m : 

9i 

92 -92 ^ („) 

9i 


(52) 




Po 


_ 2^0 


(m) 

9i 


in) 

9i 


Then, abbreviating 9 ^ = q\^^, = q'ff ', p^ = p )^’, and 


H) 


„(m) 


Ifj, = Eqs.(521 further simplify to 


In this case, we have to interpret the results in Eqs. (17 
181 for the bare parameters 


4°^ = 1 , 

9 f^ = 3-A, 
„(o) 


(0) 1 
Pi 

6°^ = 0 


= 2i-l-A (t^2), 

(*^ 1 ), 

(* ^ 1 ), 


(51) 


reflecting the fact that all sites in HNNP have a 
hierarchy-dependent, increasing degree of di = 2i — 1, 


_ c\^ 

r^i+i =r^ + 2-2 




Pn+i — V 


Pj+Pn 


^iL + l = 2 


9m' 


(53) 


Evolving the RG for k — 2 steps results in a reduced 
network that consists of 4 sites, formerly in 0, 2^“^, 2*“^, 
32 ^- 2 ^ Hence the determinant for HNNP reads 


det 


L(np)_ai 


j{k 2 ) 


det 


rk-2 —Pk-2 —24-2 —Pk-2 

—Pk-2 9fe-2 —Pk-2 0 

— 24-2 —Pk-2 rk-2 —Pk-2 

—Pk-2 0 —Pk-2 9/0-2 J 


(54) 


Jk-2) 


n -2 


{ —9/C-2 (24-2 +/'/c-2) [4Pfc_2 + 9/C-2 (24-2 —rk- 2 )] } ■ 


For HNNP, Eq.(|5^ has a fixed point at /j, — )■ oo, with tion for small A on Eq. (53) then yields [7] 
9 oo = 5, Too = 33/5,Poo = 2,/oo = 4/5. A simple perturba- 
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P>^‘ 

V ' 


26 


5 + A2^Qi 

2 -A2 ^q4 , 

5 

4 4 


(55) 


Here, ijf'^ in Eq. (521 satisfies the recursion, = 


1 

'^-2 L 

for 


j-{k 2 ) 


, which yields x ~ 2.949008159 at A: = 40 


r(fc- 2 ) 


n -2 


. Inserted into Eq. (|54|), we obtain 


(up to a factor): 
det 


L|f^)-Al 


AfV2.949008159^, (56) 


or, for the number of spanning trees in HNNP: 


n{NP) 

TTST 


' 2.949008159 


N 


(57) 


HNNP without backbone: This case is interesting be¬ 
cause HNNP without its backbone looses all its loops and 
decomposes into a fixed number (here, four) of separated 
but extensive trees [^. As Eq. ( [T7| suggests, starting 
with = 0 and = 1 for alH > 1 implies that 

= 0 and = 1 for all p and i > 1. Then, 

using the fact that = 2 throughout for i>2, 

Eq. (IT7|) collapses to 


4^+^)=4"^ ■ 


(58) 


- 2 - 


(m) 

9i 


Note that only on-site parameters renormalize, as can 
be expected for a tree. While these recursions are non¬ 
trivial, the initial conditions qf^l = 1 — A, simply lead to 

9 i ^2 “ 1 P when A = 0. Hence, so is = 1 for 

all p. 

Including 0(A)-corrections, Eq. (58 1 yields "^1 + 
Qi_ 22 ^A for large p. In the final step of the RG, these 
trees reduce to four isolated sites, at n = 0 , 2 *“^, 2 ^“^, 
and 32^“^, so that its Laplacian determinant merely has 
non-zero diagonal elements of the form q ['‘2 — 1 A^, 

where the —1 arises from the fact that these four sites 
each initially had one less link than expected from their 

level in the hierarchy. Thus, det[Lk —Al] (A^) , re¬ 
flecting the expectation that each of the four discon¬ 
nected trees merely contributes a unit factor to the count 
of spanning trees. 


( 7 |°) = 3-A 

gW=4i-3-A (i^2) 
{i^l) 

ir = 1 (* ^ 1 ) 


(59) 


( 0 ) 1 
Pi 

,(o) 


reflecting the fact that all sites in HN 6 have a hierarchy- 
dependent, increasing degree of di = 4z —3(2 ^i <k) with 
average 6 . The difference between qi and gj+i is constant 
throughout, here qi+i — qi = A. Only the renormalization 
of qi and q 2 evolve nontrivially, as before. Again, all pi 
are non-zero, encompassing the backbone links (z = 0 ) 
and all levels of long-range links (z ^ 1). But it remains 
Pi = 1 for z ^ 1 at any step p of the RG, in particular, 
p[^^ = I throughout; only the backbone po renormalizes 
nontrivially. As for HN5, in HN 6 bare links of type U are 
present in this network. Though, only li renormalizes, as 
in HN5, while li = l remains unrenormalized for all z > 2. 
Applying these considerations to Eqs. ( flTp^ results in: 


II 

9k’ 

=,k' 

9 

(m)1 

Po 


„(m) 


„(m+i) _ An) I Jl 

Pq — -I- 


(u) 

Po 


‘-kl 


in) 

Po 


in) 

9i 


(60) 


(m) 

9i 

in) 

Then, abbreviating = q['^\ = q!f\ Pfj, = Pq'\ 


and 


L = Eqs. (60) reduce to 


9m+i = 2 


9m 


rf_L+i = r^-l-4-2 


pI + 1 

9m 


(61) 


Pu+i — T 


Pn+Pu 


9m 

Z^+i = l + 2^. 

9m 


E. Case HN6 


In this case, we have to interpret the results in 
Eqs. ([I7| for the bare parameters 


For HN 6 , Eq. ( | 6 l[ ) has a fixed point at p, —>■ oo with q^o = 
5 -k 2V5, Too = 7 -k 14/V5, Poo = 2 -k V5, and Zoo = 1 + 
A simple perturbation for small A on Eq. (61) 
then yields [7] 
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g^~5 + 2V5 + A2'^Qi, 

„ 14 14-4^5 

5 ’ 


(62) 


P^~2 + 75-A2'^Qi 
~ 1H— -j= — X2^Qi 

v5 


10-3v^ 
10 ’ 
12 - 575 
To ■ 


Again, in Eq. (601 satisfies the recursion, 77^^ = 


1 

7^2 L 
for 


j(k 2) 


fe+1 

i - 

, which yields x « 4.0977251445 at fc = 40 


r{k 2) 


-2 


x^. Inserted into Eq. (541, which re¬ 


mains formally valid for HN6, we obtain (up to a factor): 


det 


Lf-A1 


A7V 4.0977251445^, (63) 

or, for the exponential proliferation of spanning trees, 

4.0977251445^. (64) 


results. Generally, the addition of extra links results 
quite naturally in an increase of complexity, measured 
in terms of the entropy-density defined in Eq. (|4]), such 
as when we progress from HN2 to HN3 and to HN5, or 
from HN2 to HNNP and to HN6. However, it is notable 
that HNNP, merely an average degree-4 network, has a 
higher complexity than HN5. Each represents a seem¬ 
ingly small variation in the basic design pattern of their 
hierarchical structure. Yet, it renders HNNP and HN6 
non-planar while HN3 and HN5 remain planar, which 
may explain the slower growth of the later with respect 
to the former. HN3 possesses the weakest growth. It 
is the only one with a small, regular degree of sites and 
with an average distance between sites that grows faster 
than logarithmic with system size (~ 7 n)- Both facts 
help suppress the combinatorial proliferation of alterna¬ 
tive paths between sites, in comparison with, say, a ran¬ 
dom regular graph of degree 3 that is locally tree-like, 
which has s = In ^ = 0.836988... [JT]. Since the num¬ 
ber of spanning trees is a metric of well-connectedness of 
networks, dynamical properties such as synchronizability 
is anticipated to be progressively better for those Hanoi 
networks of higher degree, with an additional advantage 
for those that are non-planar. 


V. CONCLUSIONS 

We have presented a RG procedure to obtain spec¬ 
tral properties of lattice Laplacians. We have applied 
this procedure to the exactly renormalizable Hanoi net¬ 
works, for which we obtain a rich set of recursion equa¬ 
tions that lend themselves to many potential applica¬ 
tions m- Here, we have analyzed these equations to 
count the asymptotic growth of spanning trees on these 
networks, and we have checked their validity for vari¬ 
ous extreme limits where they reproduce the expected 


Table I. Entropy-densities from Eq. (|^ for spanning trees on 
Hanoi networks. 


Network 

S = 7ln(#sr) 

HN2 

0 

HN3 

0.7026018588 

HN5 

1.01337412 

HNNP 

1.0814688965 

HN6 

1.4104319769 
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